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Abstract: We use the conservation law of the stress-energy and spin tensors to study the 
motion of massive zero-size objects in Riemann-Cartan geometry. The resultant world line 
equations turn out to exhibit a novel spin-curvature coupling. In particular, the spin of 
the Dirac particle does not couple to the background curvature. This is a consequence of 
its truly zero size which consistently rules out the orbital degrees of freedom. As a test of 
consistency, the wave packet solution of the free Dirac equation is considered. It is shown 
that the wave packet spin and orbital angular momentum disappear simultaneously in the 
zero-size limit. 
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1. Introduction 

The problem of particle motion in backgrounds of nontrivial geometry is usually addressed 
by using some form of the Mathisson-Papapetrou method [j], g]. One starts with the co- 
variant conservation law of the stress-energy and spin tensors of matter fields, and analyzes 
it under the assumption that matter is highly localized. In the lowest, single-pole approxi- 
mation, the moving matter is viewed as a point particle. In the pole-dipole approximation, 
its non-zero size is taken into account. 

The results found in literature can be summarized as follows. Spinless particles in the 
single-pole approximation obey the geodesic equation. In the pole-dipole approximation, 
the rotational angular momentum of the localized matter couples to spacetime curvature, 
and produces geodesic deviations ]!], ]2|, ^, |], |5[ . If the particles have spin, the curvature 
couples to the total angular momentum, and the torsion to the spin alone [y, 0, ||, ||, [l0|| . 

What we are interested in is a consistent single-pole analysis of spinning particles in 
spacetimes with curvature and torsion. This is motivated by the observation that single- 
pole approximation eliminates the influence of particle thickness, and allows the derivation 
of the pure spin-curvature coupling. In fact, this is the only way to see the influence 
of curvature on the spin part of the total angular momentum. The ambiguous algebraic 
decomposition of the total angular momentum into spin and orbital contributions are of 
no help. What we need is a truly zero-size object. As it turns out, the existing literature 
on the subject does not have this sort of prediction. 

The results that we have obtained are summarized as follows. Trajectories of spinning 
zero-size massive particles generally deviate from the geodesic lines. The deviation is due 
to the spin-curvature and spin-torsion couplings. These turn out to be different from what 
has been believed so far. In particular, the spin of the Dirac point particle does not couple 
to the curvature. If it is viewed as a wave packet solution of the Dirac equation, it does not 
couple to the torsion either. In fact, the wave packet spin and orbital angular momentum 



disappear simultaneously in the zero-size limit. We can say that Dirac point particles 
behave as spinless objects. 

The layout of the paper is as follows. In section 0, we define the conservation law of the 
stress-energy and spin tensors, and introduce the necessary geometric notions. The alge- 
braic part of the conservation equations is solved in terms of the independent variables — the 
spin tensor and the symmetric part of the stress-energy tensor. After the brief recapitula- 
tion of the covariant multipole formalism, we define the single-pole approximation for the 
independent variables, only. Section || is devoted to the derivation of the particle world 
line equations. The actual derivation is only sketched, as the method has already been 
analyzed in detail in [O]. The resulting equations of motion are compared to the pole- 
dipole equations found in literature ||, [J. As it turns out, they coincide up to a constraint 
that fixes the form of the spin tensor. This constraint is a consequence of our single-pole 
approximation, and has striking consequences on the dynamics of the Dirac particle. In 
section |||, we discuss the important case of totally antisymmetric spin tensor, and obtain a 
surprising result that spin of the Dirac particle does not couple to the background curva- 
ture. To check the consistency of our single-pole approximation, the wave packet solution 
of the free Dirac equation is analyzed. It is demonstrated that the wave packet spin and 
orbital angular momentum disappear simultaneously in the zero-size limit. In section [5], 
we give our final remarks. 

Conventions in this paper are the following. Greek indices from the middle of the 
alphabet, /jl, v, . . . , are the spacetime indices, and run over 0, 1, 2, 3. The indices from the 
beginning of the Greek alphabet, a, /?,... , take values 1, 2, 3. The spacetime coordinates are 
denoted by x^, the generic metric is denoted by g^ u (x), and r\p V stands for the Minkowski 
metric. The signature convention is ( — h +- 1-). 

2. The single-pole approximation 



We begin with the covariant conservation of the fundamental matter currents — stress- 



energy tensor t^ u , and spin tensor a x 



'p v + T\ A ) r% = t%T% v + l -o v P°K puliU , (2.1a) 

D v + T x v ^a v pff = T pff -T !Tp . (2.1b) 

Here, D u is the covariant derivative with the nonsymmetric connection r X p U , which acts on 
a vector v 11 according to the rule D^v^ = d v v^ + r^\ v v x . The torsion T x pu , and curvature 
T^vpa are defined in the standard way: 

/t-a _ 7-iA _ pX <r~ip, = ft rv _ f) r^ _i_ r^x r-^ — r^^ r-^ 

1 [iv — L Ufi 1 fiu > ' v upa — <->pi. VC j u a i V p -[- l \pi v< j l Xcr 1 up- 

The derivative D\ is assumed to satisfy the metricity condition, D\gp U = 0. As a conse- 
quence, the connection F p V is split into the Levi-Civita connection { }, and the contor- 
sion K 



pi' ■ 



1 fiu — \pu) + "■ M"' P- v — r, W ^ *v n + lp,v 



We shall also introduce the Riemannian covariant derivative V^ = D p (r — > {}), and the 
Riemannian curvature tensor R p upa = lZ p upa (r — » {}). The relation connecting the two 
curvature tensors reads: 



K^xp = R^uXp + 2V [x K^ p] + 2K\ [x K\ p] , 



where the indices in square brackets are antisymmetrized. 

Given the system of conservation equations ( |2,1[) , one finds that the second one has no 
dynamical content. Indeed, the antisymmetric part of stress-energy tensor is completely 
determined by the spin tensor. One can use ( [2.1b ) to eliminate t^ u ' from the equation 



(2.1a), and thus obtain the conservation equation, in which only r^ pv > and a Xpu components 
appear. The resulting equation reads: 

V, ( T M + X -K x fa vX P - K^ Xp a px ^ - V p a^ v A = ^a upX V p K px » . (2.2) 

This will be the starting point of the derivation of the particle world line equations. 

Let us now introduce the multipole formalism, which is necessary for the derivation. It 
has been shown in Refs. |ll|, |12J that an exponentially decreasing function can be expanded 
into a series of (^-function derivatives. For example, a scalar V(x), well localized around 
the line M, can be written in a manifestly covariant way as 



V{x) = / dr 
'M 



M(r / 4){X ^ Z) - V p (m p (t 




(2.3) 



Here, Ai is a timelike line x p = z p (t) parametrized by the proper distance, dr 2 = 
Qfiudz^dz^, and the coefficients M(t), M p (t), ... are spacetime tensors called multipole 



coefficients. It has been shown in Ref. [12| that one may truncate the series in a covariant 
way in order to approximate the description of matter. Truncation after the leading term 
is called single-pole approximation, truncation after the second term is called pole-dipole 
approximation. The physical interpretation of these approximations is the following. In 
the single-pole approximation, one assumes that the particle has no thickness, which means 
that matter is localized in a point. All higher approximations, including pole-dipole, allow 
for the nonzero thickness, and thus, for the nontrivial internal motion. 

Apart from being covariant with respect to diffeomorphisms, the series ( |2.3| ) possesses 
two extra gauge symmetries. The first is a consequence of the fact that that there are 
redundant coefficients in this decomposition. Indeed, only three out of four ^-functions in 



each term of the multipole expansion (2^3) are effective in modeling particle trajectory in 4- 
dimensional spacetime. The extra <5-function and the extra integration are introduced only 
to covariantize the expressions. The derivatives parallel to the world line are integrated out, 
as they should, considering the fact that matter is not localized in time. As a consequence, 
the parallel components of the multipole coefficients M p , M pX , ... effectively disappear. 



It has been shown in Ref. [12 1 that the corresponding gauge symmetry, named extra 



symmetry 1, in the pole-dipole approximation reads: 

5iM = Ve, 5 1 M p = u p e. 



Here, u p = dz p /dr is the particle 4-velocity, e(r) is a gauge parameter, and V stands for the 
Riemannian covariant derivative along the particle trajectory ( Vi> ^ = dv^/dr + { j^ }t> u p ). 
We see that the parallel component of M p transforms as 5\ (M p u p ) = — e , and can be 
gauged away. In fact, one can show that the parallel components of the higher multipoles 
are also pure gauge. In the gauge fixed multipole expansion, the only derivatives that 
appear are those orthogonal to the world line. 

The second extra symmetry stems from the fact that the choice of the line x^ = z^^) 
in the expansion (|2.3| ) is arbitrary. If we use another line, let us say x^ = z' p (t), the 
coefficients M, M p , ... will change to M', M' p , ... while leaving the scalar function V(x) 
invariant. The transformation law of the M-coefficients, generated by the replacement 
z p — » z' p , defines the gauge symmetry that we call extra symmetry 2. 

The extra symmetry 2 is an exact symmetry of the full expansion ( |2.3| ), but only 
approximate symmetry of the truncated series. In the pole-dipole approximation, it has 
the form 

5 2Z p = e^ , 6 2 M = Mu p Ve p , S 2 M P = -Me p , 

provided the M-coefficients are subject to the hierarchy 

M = O , M p = £>! , M px = 2 , 

and the free parameters e p (r) satisfy e p = 0\. Here, O n stands for the order of smallness, 
and the condition e M = 0\ ensures that the order of truncation is not violated by the 
action of the symmetry transformations WA - I n the pole-dipole and higher approximations, 
fixing the gauge of extra symmetry 2 defines the particle centre of mass. In the single-pole 
approximation, the extra symmetry 2 is trivial. 

Now, we shall replace the general function V{x) with the stress-energy and spin tensors 
of the localized matter. In order to describe a strict point particle, we choose T^ pv > and 
g pv in the form 

t (a-) = f dTb p »(T) 6(D) (*Z Z \ (2.4a) 

Jm v^g 



a^= I dTC Apv {T)- ^^^, (2.4b) 

Jm V-9 

where 6^"(r) and c Xiiu (t) are the corresponding multipole coefficients. We emphasize here 
that this is not how single-pole approximation is defined in the existing literature [||. |J. 
There, the antisymmetric part of stress-energy tensor r pu has also been treated in the 
single-pole manner. As T"*"' is not an independent variable, this imposed unnecessary con- 
straints on a xpu . In particular, the spin of the Dirac particle was ruled out. To overcome 
this problem, the authors of Ref. || abandoned single-pole in favour of pole-dipole ap- 
proximation. Their subsequent limit of vanishing orbital angular momentum should have 
brought them back to the single-pole regime. In what follows, however, we shall demon- 
strate that it is not quite so, and that such a limit is not equivalent to the single-pole 



approximation as defined in (2.4) 



3. Equations of motion 

The particle equations of motion are derived in the following way. We insert ( p.4| ) into 
Q2.2p , and solve for the unknown variables z(t), b pu {r) and c Xiiu (t). The algorithm for 
solving this type of equation is discussed in detail in pT| , |l2|, and here we only sketch it. 
The first step is to multiply the equation (2.2) with an arbitrary spacetime function fn(x), 



and integrate over the spacetime. The resulting equation depends on the function /„ and 
its first and second covariant derivatives, evaluated on the line x p = z p (t): 



Jdr[c^"f^ p + (lT - K^ Xp c^ + l -K X/ rc^)f,, v + \c vpX {v»KP Xv )f, 



0. 



where f K „ = {S/ V fp)x=z, fp.\v P = (V p V !/ / M ) a=z . Owing to the arbitrariness of the function 
fp(x), the terms proportional to its independent derivatives separately vanish. To find 
the independent derivatives of the test function /„, we make use of the particle 4- velocity 
u p = dz^/dr, and the Riemannian covariant derivative along the particle trajectory V. 
The 4- velocity u p is normalized as u p u p = — 1, and the action of V on a vector field v p (t) 
is defined by Vv^ = dv^/dr + {j^ }u u p . Next, we decompose the derivatives of the vector 
field ffj,(x) into components orthogonal and parallel to the world line x^ = z^(t): 

U;X = fp\ ~ ^aV/ m , (3.1a) 

f K (Xp) = fpXp ~ 2h p(\ u p) + h pU X u p , (3.1b) 

fp;[Xp] = l^R a pXpfv ■ (3.1c) 

Here, the orthogonal and parallel components are obtained by using the projectors 

P^ = Si + u p u u , PI = -uPuv. (3.2) 

More precisely, f^ x = P± x f^, f pXp = P ± lP ± u p f^ au) , h^ x = F \ a x u v 'j ' p . {au) and h p = 
u a u v ' f^cju)- Direct calculation yields 



K =w/ M -(v«")/^, 

K P = P^fpu ~ (Vn p )V^ + \p4u v R% vX U , 

which tells us that the only independent derivatives on the line x p = z p (t) are f p , f pu and 
fjiup- We can now use (34.) and ( |3.3D to group the coefficients in terms proportional to the 
independent derivatives of / M . The obtained equation has the following general structure: 



dr 
M 



^7i + ^"7i + n 



o. 



where X^ up , X pv and X p are composed of various combinations of multipole coefficients 
b^ and c XfJ,1/ , external fields K x pu and R p vp(J , and their derivatives. In all the expressions, 
the external fields are evaluated on the world line x p = z p {t). Owing to the fact that 



ffii fuu an d fuvp are independent functions on the world line, we deduce that the X-terms 
must separately vanish. The equation X fiup = has a simple algebraic form 

which is easily solved for c pu . This yields 

where s pl/ = —s u>1 and s pv = —s^ = s u p are totally antisymmetric, but othervise 
free parameters. The equations X pv = and X^ = are much more complicated. The 
procedure goes as follows. First, we use the above decomposition of c xpu to perform a 
similar split of the b pv coefficients. A new free parameter m{f) appears to characterize 
the leading term of b pu . Then, the equations X^ u = and X p = are rewritten in terms 
of the undetermined parameters m, s pu and s '"', and properly rearranged. Skipping the 
details of the diagonalization procedure, which has thoroughly been demonstrated in Ref. 



12], we display the final result: 
• the world line equation 

V 



m/ + 2u p (Vs pp + D pp ) ] - u v s pa R» vpa = ^c vpX V p K pXv , (3.4a) 



• the spin precession equation 

P^pPZ (Vs pa + D pa ) = 0, (3.4b) 

• the stress-energy coefficients 

IT = mu p u u + 2u x u ( - fl Vs^ x - -K Xp ( p c u ^ xp , (3.4c) 

• the spin tensor coefficients 

c Xpv = 2u x S p u + s Xpu . (3.4d) 

In these equations, the scalar m(r), and the totally antisymmetric tensors s pu (t) and 
s pu (r) are free parameters. They determine the stress-energy and spin tensors via fl3,4dj 
and fl3.4d|) . The shorthand notation 

D ^u s K [n XpC P^} + \k x /c u]pX 

is introduced to simplify the cumbersome expressions. 

The obtained single-pole equations differ from the known pole-dipole result || [7|, |8|, 
pi, |10| by the presence of the constraint (3.4d|). It is a consequence of our assumption 



that the particle has no thickness, and therefore no orbital degrees of freedom. In the 
existing literature, an analogous but more restrictive constraint appears in this regime 
0, J9|. This is because the antisymmetric part of the stress-energy tensor T"*"! has been 
treated as an independent variable, in spite of the restriction ( p. lb ). In our approach, the 



only independent variables are a pv and r^ v \ and the resulting constraint ( 3.4d| ) is not 



so strong. In particular, it does not rule out the free Dirac field, or any other massive 
elementary field. Indeed, the formula c XfJ,u = u x s fJ-u -\ — u^s u '^ for the spin tensor of the 
elementary particle of spin s (see Ref. [10|) is a special case of ( |3.4d ). 



In what follows, we shall examine the special case of spin 1/2 pointlike matter. Surpris- 
ingly, we shall discover that spin 1/2 does not couple to the curvature, leading to geodesic 
trajectories in torsionless spacetimes. 

4. The Dirac particle 

The basic property of Dirac matter is the total antisymmetry of its spin tensor a^ v . As a 
consequence, the coefficients c A ^ are also totally antisymmetric, and the constraint ( |3.4d ) 
implies 

a"" = 0. (4.1) 

The vanishing of the s^ v component of the spin tensor has far reaching consequences. 
First, we see that the spin-curvature and spin-orbit couplings disappear from the world 
line equation ( 3.4a| ). Second, the spin precession equation ( 3.4b ) becomes a constraint 



equation. If we define the spin vector s^ by s^ vp = e fJ,up s\, and the axial component of 
the contorsion K p as K^ = e pvi K up \, where e^t" 1 is the covariant totally antisymmetric 
Levi-Civita tensor, the equations (3.4) become 



V (mu» + K^s^u v \ + }-s p V^K v = 0, (4.2a) 

K^s v l = 0. (4.2b) 

As we can see, the spin couples only to the axial component of the contorsion, which 
means that Dirac point particles follow geodesic trajectories in torsionless spacetimes. At 
the same time, the absence of torsion trivializes the equation (4.2b), and no information on 



the behavior of the spin vector is available. If the background torsion has nontrivial axial 
component, a geodesic deviation appears, but also a very strong constraint on the spin 
vector. Indeed, the equation ( [4.2b| ) implies that the orthogonal component of s^ always 
orients itself along the background direction K*?. This unusual behavior suggests that the 
spin vector of the Dirac point particle might be zero, after all. In fact, the world line 



equations (|3.4j ) are derived under very general assumptions of the existence of pointlike 
solutions in an arbitrary field theory. They do not care about peculiarities of specific 
theories or specific types of localized solutions. In what follows, we shall analyze the wave 
packet solutions of the flat space Dirac equation, with the idea to check if they can be 
viewed as point particles. 

Let us construct an example. We start with the free Dirac Lagrangian 

C = % - [fifdrf - (fy^/V)] - rmfjiP , 

where Dirac 7-matrices satisfy the usual anticommutation relations {7^,7^} = —2rj fJ ' h ', and 
are used in their conventional representation (75 = i7°7 1 7 2 7 3 ). Then, we construct a wave 



packet. The wave packet is a solution which is well localized in space, but resembles a plane 
wave inside. To be viewed as a particle, its size £ is considered in the limit £ — > 0. At the 
same time, the particle stability is achieved in the limit X/£ — > 0, where A is its wavelength. 
We construct it as follows. At the initial moment x° = 0, we choose the configuration 

^(r,0) = Ae &ip p (r,0), 



where 



i/j p (x) 



k° + 



in 



2m 



1 


fc 3 

W+m 





Ak u x^ 



is the plane-wave solution of the Dirac equation (i^^d^ — m) tp = 0. It propagates along the 
x 3 -axis {k l = k 2 = 0, k° = ^Jm 2 + (/c 3 ) 2 ), and is polarized upwards, for convenience. The 
exponential function exp(— r 2 /£ 2 ) cuts out a small piece of the plane wave, and defines its 
size £, while A is the overall amplitude of the packet. The wavelength A is proportional to 
1/|A;|. Using the Dirac equation, we can calculate time derivatives and thereby determine 
time evolution of this packet. In fact, we only need first time derivatives, as neither t^ u 
nor a ^ v depend on higher derivatives: 



I jJLU 



i [V^aA^ - {d„ip) ltii>] - ZrfuvC . 



a 



Xfii/ 



T X/iup 



^757 P V • 



(4.3) 



The wave packet expressions of these currents at x° = are obtained straightforwardly: 



-(00) 



-2\A\ z e~ 



m 



_(0«) 



-2\A\ 2 e~ 



k° 



hi 



k 3 T] 3a 



,k 3 



a ™=i A \ 2 e~ir-, a 012 



-2\A\ 2 e 



X P a/33 

£ 2 



-\A\ 2 e~ 



in 



,k° 



(4.4a) 
(4.4b) 

(4.4c) 



m m 

where only non-vanishing components are displayed. Now, we want to rewrite the currents 
( |4.4| ) as a series of 5- function derivatives. We first fix diffeomorphisms by imposing the 
condition g^ u = rjavi an d extra symmetry 1 by keeping only spatial components of the M 
coefficients (M° = M op = • • • = 0). The decomposition formula fl2.3| ) is thereby reduced to 

V(x) = I dr \m5^ (x-z)- d a (M a S^ (x - z)) + ■ ■ ■ 

In general, the line x M = z^(j) is arbitrary, but the simplest expressions are obtained if it 
coincides with the wave packet trajectory. Thus, we choose z 1 = z 2 = in accordance with 
the the fact that the packet propagates along the x 3 axis. As for the z 3 component, we do 
not need its full r dependence because we are only interested in the packet behaviour at 
x° = 0. There, the proper length t is chosen in accordance with z°(0) = z 3 (0) = 0, which 
is sufficient for the proper definition of the (5-expansion at x° = 0. 



The multipole coefficients are obtained by multiplying V(x) by a number of (x a — z a ) 
factors, and integrating over the 3-space. In our example, the simple integration of (f4.4| ) 
yields the monopole coefficients b^ v and c Xllu , while multiplication with (x a — z a ), and 
subsequent integration gives the dipole coefficients W va . The resultant non-zero monopoles 
are 

b 00 = a£ 3 (k ) 2 , b 33 = a£ 3 (k 3 ) 2 , b 03 = a£ 3 k°k 3 , 

C 123 = Vlfe 3 , C 012 = -Pk°, 

2 2 

while there are only two non-zero dipoles, 

6 012 = _ 6 021 = ^3 fc 0_ 

4 
The higher multipoles are of the order £ 5 (fc 3 ) 2 or higher. Here, a(k) is the overall factor 
whose explicit form is not needed in the subsequent discussion. 

Now, we shall consider the single-pole limit £ — > 0, while respecting the wave packet 
stability condition A <C £. First, we choose the overall amplitude a(k) in the form 

A 2 

thereby normalizing the monopole coefficients to be of the order of unity. Then, using the 
single-pole behavior k 3 ~ k° ~ A" 1 , we find 

6°o „ b 33 ~ 6 03 ~ 1 , c 123 ~ c 012 ~ A , 6 012 ~ A . 

The higher multipoles are of the order £ 2 , and thus, neglected. This is a realization of the 
pole-dipole approximation. In the single-pole regime, however, only the lowest terms are 
retained in the limit £ — > 0. This means that, respecting A <C £, the terms proportional to A 
must also be dropped. As a result, the spin monopole coefficients c '*" vanish simultaneously 
with the orbital dipole coefficients b^ . 

The reason for this unusual behaviour is found in the constraint c 012 = 2 ft 012 . It is 
obtained by the integration of the more general relation 



that is found to constrain the wave packet currents (4.4). It relates the wave packet spin 
to its orbital angular momentum, so that the expected disappearance of orbital degrees of 
freedom in the limit £ — > is followed by the unexpected disappearance of the spin itself. 
To summarize, we see that the spin vector s^ vanishes in the single-pole approximation, 
and the particle trajectory becomes a geodesic line even in the presence of torsion. The 
validity of this conclusion, however, demands some sort of equivalence principle to hold. 
This is because the considered wave packet is a solution of the free Dirac equation, and 
the inclusion of curvature or torsion may destroy it. What we can do is to consider weak 
gravity, so that terms quadratic in curvature and torsion are neglected. In that case, the 
free wave packets are a good approximation to the exact solution, which implies that Dirac 
point particles behave as spinless objects in an external gravitational field. They can still 
probe the spacetime curvature, but for the probe of the background torsion, one needs a 
thick particle. 



5. Concluding remarks 

In this paper, we have considered the motion of point particles with nonzero spin in space- 
times with curvature and torsion. Using the covariant multipole formalism developed in 



Ref. [12], the world line equations are derived in the lowest, single-pole approximation. 
This way, the particle thickness, and the corresponding internal motion, have been elimi- 
nated. Only mass and spin remained to characterize the particle internal structure. 

In our approach, the single-pole behaviour has been adopted for truly independent 
variables only. In particular, the antisymmetric part of the stress-energy tensor has been 
eliminated from the conservation equations, prior to imposing the single-pole regime. As 
a consequence, our single-pole analysis turned out to differ from the existing literature. 

The obtained equations of motion are found to differ from the known pole-dipole 
equations by the presence of a novel constraint on the particle spin tensor. With this 
constraint, the spin of the Dirac point particle turned out not to couple to the background 
curvature, leading to geodesic trajectories in torsionless spacetimes. In the presence of 
torsion, however, a geodesic deviation appears, but also a strong constraint, suggesting 
that the spin of the Dirac point particle might be zero. Being a consequence of our single- 
pole approximation, this unusual result has been checked by the explicit construction of 
the zero-size Dirac particle. To this purpose, a wave packet solution of the Dirac equation 
is considered in the limit of small size and wavelength. The expected single-pole behaviour 
has been verified, but also, the spin tensor has been found to disappear in this limit. In 
an attempt to explain this unusual behaviour, a relation between spin and orbital angular 
momentum has been discovered to hold in our wave packet example. 

Before we close our exposition, let us comment on the possibility that the disappearance 
of spin in the zero-size limit might be a general property of all point particles. First, we 
notice that there is one more conserved current in the Dirac theory — the U(l) current 
j^ = tpj^i/j. In our wave packet example, it is proportional to the wave vector fc M , and is 
related to the stress-energy t^ v . By a close inspection, we find that the following manifestly 
covariant relation holds: 

T^ oc j", (5.1) 

Its physical meaning is best seen in the rest frame where it reduces to r Q ° = j a = 0. It 
tells us that the two currents are mutually proportional, i.e. that the energy and charge 
flow in the same direction. In the limit \/i — » 0, I — > it implies the constraint (fl.5|), and 
thus, explains why Dirac point particles have no spin. 



The relations analogous to (5.1) might exist quite generally. There is nothing special 
about the statement that all the particle charges flow in the same direction. This is 
something one would expect to hold for any type of point-like matter. However, we need 
the relation ( |5.1| ) to hold for thick particles, as well. Only then, and only for Dirac matter, 
the consequence ( [4.5D has been derived. If ( [5.1] ) were a general property of the localized 



matter, the disappearance of spin in the zero-size limit might be a feature of all massive 



point particles. Indeed, when applied to other spins, the relation (5.1) implies 

x * T m _ x /3 r (oa) = a i0oa] + x w daa om + div _ 



10 



After the integration, the divergence term vanishes, while l.h.s. and the second term on the 
r.h.s. give dipole coefficients of the order 0\. In the single-pole regime, both disappear, 
and we end up with the vanishing spin. 



References 



[1] 
[2] 
[3] 
[4] 



M. Mathisson, \Acta Phys. Polon. 6 (1937) 163 . 
A. Papapetrou, \Proc. R. Soc. A 209 (1951) 248j . 



W. Tulczyjew, Acta Phys. Polon. 18 (1959) 393. 



A. H. Taub, \J. Math. Phys. 5 (1964) 112| . 



G. Dixon, Nuovo Cim. 34 (1964) 317; Nuovo Cim. 38 (1965) 1616; Proc. R. Soc. A 314 



(1970) 49S|; \Proc. R. Soc. A 319 (1970) 509| ; \Gen. Ret Gray. 4 (1973) 199fc \Proc. R. Soc. A 
319 (1974) 509J . 



[6 

[T 

[8] 

[9] 

[10 

[11 
[12 



A. Trautman, 



Bull. Acad. Pol. Sci., math, astr. phys. 20 (1972) 895 



F. W. Hchl, Phys. Lett. A 48 (1976) 393 



P. Yasskin and W. Stoeger, Phys. Rev. D 21 (1980) 2081 



K. Nomura, T. Shirafuji, and K. Hayashi, \Prog. Theor. Phys. 86 (1991) 1239| . 
K. Nomura, T. Shirafuji, and K. Hayashi, \Prog. Theor. Phys. 87 (1992) 1275| . 



M. Vasilic and M. Vojinovic, Phys. Rev. D 73 (2006) 124013, gr-qc/0610014 



M. Vasilic and M. Vojinovic, J. High Energy Phys. 07 (2007) 028, gr-qc/0707.3395. 



11 



